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The conductance of an open quench-disordered two-dimensional (2D) electron system subject to 
£\j , an in-plane magnetic field is calculated within the framework of conventional Fermi liquid theory 

applied to actually a three-dimensional system of spinless electrons confined to a highly anisotropic 
(planar) near-surface potential well. Using the calculation method suggested in this paper, the mag- 
netic field piercing a finite range of infinitely long system of carriers is treated as introducing the 
■ additional highly non-local scatterer which separates the circuit thus modelled into three parts — 

£f~^ ' the system as such and two perfect leads. The transverse quantization spectrum of the inner part 

of the electron waveguide thus constructed can be effectively tuned by means of the magnetic field, 
even though the least transverse dimension of the waveguide is small compared to the magnetic 
length. The initially finite (metallic) value of the conductance, which is attributed to the existence 
of extended modes of the transverse quantization, decreases rapidly as the magnetic field grows. 
. This decrease is due to the mode number reduction effect produced by the magnetic field. The clos- 

r/) ing of the last current-carrying mode, which is slightly sensitive to the disorder level, is suggested as 

D ' the origin of the magnetic-field-driven metal-to-insulator transition widely observed in 2D systems. 
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I. INTRODUCTION 
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The conduction properties of low-dimensional electron and hole systems with the disorder of different origin have 
long been the subject of active research. Investigations into such objects of mesoscopic size have currently become 
I ' particularly intensive in view of their applied importance (semiconductor heterostructures, quantum dot devices, etc.), 
, on the one hand, and due to the intriguing uncommonness of the obtained results, on the other. The most puzzling 
phenomenon which has not as yet been understood in full measure is the transition of 2D electron systems from con- 
ducting to insulating state (the extensive bibliography on the subject can be found in Refs. 0,0)- Observation of this 
phenomenon evidently contradicts the common view stemming from the well-known scaling theory of localization [j| . 

A diversity of ideas have been put forward to explain the "anomalous" metallic behaviour of 2D systems at extremely 
low temperatures. Among physical mechanisms assumed to be responsible for such behaviour the most appropriate is 
the Coulomb (e-e) interaction of carriers. The assumption is mainly based on the fact that the dimensionless parameter 
r s (the ratio of interaction energy to Fermi energy of electrons) which characterizes this interaction amounts to the 
value of > 10 in planar systems of Si MOSFET-type as well as in GaAs/AlGaAs heterostructures, which would suffice 
to raise serious doubts about the adequacy of Fermi-liquid approach with regard to such systems. 

Note that in condensed matter physics an important paradigm is the quasi-particle notion. A number of phenomena 
can be understood in terms of weakly interacting quasi-particles, although Coulomb interaction between the electrons 
is rather strong. The Fermi liquid description of quasi-particles has for a long time been considered to be efficient 
in low dimensions in the presence of disorder 0, Q until recent experiments have challenged this view berious 
problems have also arisen when interpreting numerous experimental data on decoherence time saturation at low 
temperature, reviewed in Ref. 0, as well as explaning unexpectedly large persistent currents in normal metals 0- 

At this stage, the lack of a comprehensive theory for strongly correlated electrons does not allow for certain 
conclusions to be made about the role of e-e interaction in the detected phenomena. Since in different existing 
theories this interaction is evaluated quite differently, viz. from promoting localization [ij llOf to preventing its 
appearance 0, [T3. IT^ | , it would be quite reasonable at first to describe the observed effects in frames of one-particle 
approach and then to regard Coulomb interaction as the added dephasing factor. Only should one fail in obtaining 
the proper result on this canonical way may it be rational to think of Coulomb correlations as a prevalent cause of 
forming the continuous component in the electron spectrum of 2D systems. 
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Previously , the one-particle theory capable of explaining the conducting ground state of weakly disordered finite- 
size 2D systems was developed in terms of electron states pertinent to open quantum wells of waveguide configuration. 
According to this theory, the conductance must exhibit non-exponential coordinate dependence at any length of the 
system provided that there exists more than one extended mode of transverse quantization (so called waveguide 
mode), commonly referred to as the open conducting channel |l5j . In Ref. |l4j it was proven that scattering by 
static random potential, even though it is of elastic nature from the viewpoint of one-electron dynamics, may well 
result in dephasing the initially coherent extended collective modes. This leads to the reduction in the primordially 
finite-value ballistic conductance rather than it grows from the "localized" exponentially small value. The primordial 
conductance is quantized in steps, each being equal to the conductance quantum Go = e 2 /-Kh, regardless of the length 
of the bounded 2D system of carriers. Inclusion of imperfections allowing for scattering between different extended 
waveguide modes results in the conductance decrease from ballistic to diffusive value. The expression for the diffusive 
conductance tends asymptotically to the standard Drude form if the quantum waveguide is sufficiently wide. The 
above scenario of the conducting state of a real 2D electron system is evidently opposite to the one stemming from 
well-known scaling theory, in which the metallic behaviour of the conductance can only result from some kind of 
dephasing of the electron states initially localized by disorder. 

Using the approach taken in Ref. |b|, in Ref. [lfj one-particle theory of the metal- insulator transition (MIT) 
in 2D systems was developed. This theory provides a deeper insight into the peculiarities of the effect observed 
experimentally. According to 0], (i) MIT detected in Si-MOSFETs as well as in structures of GaAs/AlGaAs-type 
should be regarded as a true quantum phase transition which takes place both in disordered and perfect systems; 
(ii) straight from the metallic side of the transition the conductance reaches the value equal, by the order of magnitude, 
to the standard conductance quantum (theoretically, in the limit of an ideally perfect system, the conductance jumps 
exactly by Go as the last conducting channel is closed). 

The results [bH 1 1 6| were obtained in the Fermi-liquid description of the electron system, without taking explicitly 
into account Coulomb interaction of carriers as well as spin effects. Meanwhile, the very fact that the extended 
waveguide modes were used as the initial quasi-particles (which have unlimited spatial extent, in contrast to Anderson- 
localized states) permits hoping that the results of Refs. will hold true, at least qualitatively, even on condition 
that Coulomb correlations are also taken into account. 

Although the model of 2D MIT suggested in Ref. 0] correlate well with the main experimental facts, it does not 
seem to be rather convincing to make far-reaching conclusions. The problem of great concern, which has not been 
touched upon in |16| and has not been so far explained theoretically is the unusual nature of the strong localizing 
effect produced by the relatively weak in-plane magnetic field. Normally, the magnetic field is known to suppress 
the conductance of bulk metallic samples due to coupling to the orbital motion of electrons. However, if the field is 
applied parallel to a confined system whose actual thickness is small as compared to the magnetic length, it is natural 
to expect the orbital effect to be highly suppressed. Only corrections from e-e interaction of carriers are retained in 
the case of spinless electrons Q . 

With this prevailing point of view, it came as a surprise when dramatic suppression of conductivity was observed 
in Si-MOSFETs subjected to the in-plane magnetic field flTl . The effect of this field was found to be so significant 
that it causes the zero- field 2D metal to become an insulator |l8l Il9l I20I l2l| . Inasmuch as magnetic fields appropriate 
for such a transition are strong enough to completely polarize electron spins, this prompted many of the researchers 
(see, e. g., Refs. [^HHEl) to assert categorically that it is exactly the spin polarization that should be regarded 
as a physical origin of 2D MIT in a parallel magnetic field. This point of view is also supported by the lack of a 
comprehensive theory of the electron transport in real quasi-two-dimensional (Q2D) systems, i. e. planar systems 
of finite, though rather small, thickness. Although in some papers (e. g., [2M l2*q1 the orbital effect of an in-plane 
magnetic field was, in a way, analyzed as well, the results are still not rated as fully convincing. This is probably 
concerned, on the one hand, with the fact that the relatively simple model proposed in |25j is incapable of explaining 
a quite abrupt transition over the magnetic field between metallic and dielectric regimes. On the other hand, the 
results obtained in Ref. 26| only exhibit the trend in behaviour of the Q2D system conductance with a growth in the 
magnetic field, since the applied calculation method enables one to analyze only weak-localization corrections in the 
metallic phase, thus being incapable of capturing the coarse effect such as MIT. 

In this paper, the orbital effect of in-plane magnetic field upon open quasi-2D electron systems of waveguide 
configuration is investigated using the theoretic-field method previously employed in Refs. [T3. IT^|. The advantage 
of the approach is that it enables one to take into account in a similar fashion all the one-electron states, without 
distinguishing them by the kind of trajectories, which in Q2D case can be either of sliding type (Levi flights) or those 
frequently scattered at side boundaries of the system. Although the method is widely applicable, here we will examine 
thoroughly the specific case of weak electron-impurity scattering and weak magnetic fields (the appropriate criteria 
are listed below). 

The notion of weak magnetic field implies the certain relationship between geometric parameters of the electron 
waveguide and the classical cyclotron radius of electron trajectories. At the same time, in general it does not imply low 
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mixing of the electron states. The magnetic-field-induced mode entanglement may be rather strong, but it does not 
give rise to noticeable dephasing of the transverse-quantized states unless there exists some random short-correlated 
potential as well. The most important effect of the magnetic field is that it influences significantly the mode content of 
the electron spectrum, thereby resul ting in noticeable change in the number of current-carrying modes and in spectral 
width of their energy levels. In Ref. |2j|, it was proven that the number of conducting channels becomes smaller when 
the magnetic field increases, being simultaneously accompanied by the gain in coherence of the electron states. It is 
just the reduction in the number of extended modes that can be viewed as the physical origin of the transition of 
Q2D electron system from conducting to insulating state. 

To justify the sensitivity of the mode state spectrum as the major cause of the observed MIT it is necessary to 
trace the spectrum peculiarities in terms of the observable quantity, say, the conductance. To accomplish this task, 
it would be inconvenient to use the Landauer approach, since it does not look into the coherence properties of the 
internal spectrum of the system in question, operating it as a scattering object integrally. Therefore, in this study 
we apply the linear response theory |28j which permits us to calculate the conductance on a microscopic level. 



II. THE MODEL 



Two-dimensional electron and hole systems in use, in view of their open property in the direction of current and the 
resemblance of master equation to that of the classical wave theory, can be simulated as planar quantum waveguides 
whose transverse design is governed by the lateral confinement potentials. Although near-surface potential wells in 
Si MOSFETs and GaAs/AlGaAs heterostructures are close in shape to triangular or parabolic form [l2l |29|. this 
fact is of minor importance for its principal application which is to restrict electron transport in the direction normal 
to heterophase areas, thus resulting in transverse quantization of the electron spectrum. With this consideration in 
mind, we assume, to simplify calculations, the Q2D system of carriers to have the form of three-dimensional planar 
"electron waveguide" of rectangular cross-section (see Fig. , which occupies the coordinate region 

xe (-L/2.L/2) , 

ye [-W/2,W/2] , (1) 
z e [-H/2,H/2] . 

The length L, the width W, and the height H of the waveguide will be regarded as arbitrary, within the restrictions 
imposed below. 




L 



FIG. 1: Configuration of the electron waveguide and the magnetic field. 



In practice the change of a quantum waveguide thickness implies alteration of the width of the near-surface potential 
well and, as a consequence, of sheet density of the carriers. Inasmuch as this density is known to follow the variation 
in depletion voltage under the simple law |3fl IsH ] , the results obtained below as a function of a waveguide thickness 
can be easily related to the experiment. 

We will examine the T = magnetoresistance of a Q2D electron system by expressing the dimensionless (in units 
of e 2 /irh) conductance g xx in terms of one-particle propagators. Assuming the system of units with h = 2m = I (m is 
the electron effective mass), the static conductance is given by 

[G £ Q F (r,r')-G^(r,r')] 
[G a ep (r',r)-G r eF (r',r)] , (2) 

where Ge^(r,r') is the retarded (advanced) Green function of the electrons of Fermi energy Ep, A^r) is the i-th 
component of the external vector potential A(r). Integration in is performed over the region occupied by 



g xx (L,B) = 
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the quantum waveguide, spin degeneracy is taken into account by the factor of 2. Note that throughout this paper 
the Fermi energy (or the chemical potential, as the situation requires) will be considered to have a constant value, 
independently of the confinement potential. This is undoubtedly true on the metallic side of the MIT discussed below. 
Moreover, on just-dielectric side of the transition this assertion is also valid since the "expanded" electron system, 
which includes the attached leads, in the static limit is on a homogeneous equilibrium state provided that its segment 
of interest is open, at least, partially. 

Within the model of isotropic Fermi liquid, the retarded Green function of electrons subjected to a static magnetic 
field obeys the equation (all indices are omitted for brevity) 



k% + iO- V(r) 



G(r,r') = <S(r-r') 



(3) 



where $o — hc/e is the magnetic flux quantum, k F = Sf, V(r) is the random potential due to impurities or the 
roughness of the confining well boundaries. We will specify this potential by zero mean value, ^V(r)) = 0, and the 
binary correlation function (V(r)V(r')} = QW(r — r'). The angular brackets are used for configurational averaging, 
W(r) is the function normalized to unity and falling off its maximal value at r = over the characteristic length r c 
(the correlation radius) . 

Equation © must be supplemented with the appropriate boundary conditions (BC). We will regard the electrons to 
be confined by infinitely high potential walls at side boundaries of the region and specify this fact by the Dirichlet 
conditions as follows, 



G(i 



y = ±W/2 
z = ±H/2 



(4) 



As far as open ends of the system are concerned, the BC problem is resolved somewhat less trivially. The mere 
fact that the system is open, even partially, implies non-Hermicity of the operator © in the domain QJ. This may 
cause some vagueness regarding the applicability of the formula (0) , whose derivation relies essentially on Hcrmitian 
property of the Hamilton operator. In the case of a finite-length system, the hermitizing BC at x = ±L/2 would in 
fact correspond to its closeness (or periodicity) in ^-direction, which does not conform with the requirement of current 
overflow between independent reservoirs. 

In this study, in view of the chosen Green function formalism, to specify the openness of the quantum system 
we employ the method based on the analogy between the problem @ and that of the monochromatic point source 
radiation in a classical waveguide. When solving the latter problem, Sommerfeld radiation conditions are normally 
used [32l l33l | , which imply for the source positioned at some finite coordinate the existence of solely outgoing waves at 
infinity. In order to adapt these conditions to the system under consideration it is necessary to prolong the disordered 
and magnetic-field-subjected segment Q of the electron waveguide with semi-infinite ideal leads, in which the electron 
waves generated at some point inside the segment could propagate freely to infinity, not being subjected to any kind 
of backscattering. Then, joining of the solutions to Eq. © at interfaces x = ±L/2 within the electron waveguide 
results in the complete solution corresponding to the infinite open system, thereby giving rise to the correct BC at 
the ends of the segment of interest. This somewhat troublesome procedure will be done in the Appendix, where the 
trial Green functions, which serve as a basis for the exact solution to equation are obtained. 

Apart from the openness BC, one more problem is to be resolved before we proceed to the conductance calculation. 
Specifically, this is the choice of the vector potential gage corresponding to the in-plane magnetic field B = (B x , B y , 0). 
Here, the main motivation is to avoid quadratic confinement in x-direction. The lack of such a confinement enables one 
to employ the transverse eigenfunction basis consisting of waveguide modes, both extended and evanescent, instead 
of normally used isotropic plane- wave basis. The choice of the latter basis would make the BC at side boundaries not 
easy to be satisfied, which would result in a substantially complicated solution. We prefer the gage 



A(r) = {B y z,-B w z,0) , 



(5) 



which yields equation © of the form 



G(r,r') = <J(r-r') . 



(6) 



Eq. 10 is suitable for reducing the initially three-dimensional stochastic problem to a set of more readily solved 
one-dimensional problems. Such a reduction is particularly helpful in view of the openness BC, which is most easily 
formulated in one spatial dimension. 
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III. REDUCTION TO ONE-DIMENSIONAL DYNAMIC PROBLEMS 



Clearly, the initially three-dimensional dynamic problem might be reduced to one-dimensional problems if one 
succeeded, yet hypothetically, in finding the actual set of decoupled transverse eigenmodes of the random Hamiltonian 
in Eq. @. As a matter of fact, one can do so, though indirectly, in terms of standard Fourier transformation of this 
equation over transverse vector coordinate r_L = (y, z). In our confinement model, the appropriate eigenfunctions are 



\r±;fj) = 



VWH 
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(7) 



with fi = (n,m) being the vector mode index conjugate to (n, m £ N). With functions (7J, equation © is 
transformed to the set of one-coordinate, yet strongly coupled, equations for the mode Fourier components of the 
function G(r,r'), viz. 



Gfj.fj.'ix^x') - ^ tinv{ x )G vl j,>(x,x') = S^Six - x') 



(8) 



Here 



is the unperturbed mode energy, 



/7rn\ 2 /7rm\ 2 



V W J 



H J 



(9) 



H 2 
121% 



(10) 



is the diagonal-in- mode- indices matrix element of the total potential which includes the disorder- induced part, V(r), 
and all the magnetic-field-related terms in l.h.s. of Eq. ©, Ib — \J Qq/2ttB is the total magnetic length. The term 
V/j.^x) in Eq. lfTt]|) is the diagonal element of the mode matrix ||V^u/||, whose elements are evaluated as 



V liU (x)= / dTx(r±;ti\V(r)\rx;f) , 
Js 



(11) 



integration is over cross-section 5* of the quantum well l(Tjl. Off-diagonal matrix elements U^uix) in Eq. JSJl likewise 
include the disorder- and the magnetic-field-originated potentials, 



Uy, v {x) = V^(x) + 2iH 



' q( x ) 



l 2 v dx 



C * 

(0 



(12) 



Here, partial magnetic lengths k (with i = y,z) are given by I 2 — $o/27r|Sj|, S)ll and are the model-specific 
numerical coefficients, each of the order unity, 
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-(to - mi) 
-(m — toi) 



(13a) 
(13b) 
(13c) 



In Eqs. I|13|) . Si = Bi/\B%\ = sgn£?i and the notations for mode indices are employed such that /x = (n, m), 
1/ = (ni,TOi). Note that the specific form of coefficients (II 31) does not significantly affect the algebraic structure 
of Eq. © and hence the subsequent transformations. This fact is the main cause for the final results being uncriti- 
cally sensitive to the actual form of the confinement potential. 

In Eq. JS}, matrix elements Vfi(x) and U^ix) may be regarded as the effective potentials responsible for intra- mode 
and inter-mode scattering, respectively. We thus adopt the approach where interactions of the electron system with 
both the disorder potential and the magnetic field are exploited similarly, i. e. they are described by the additive 
static potentials which are basically different in their correlation properties. 
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To proceed further with the transverse mode decoupling, we will follow the procedure outlined in Refs. 0,0]. It 
was shown there that the off-diagonal components of the mode Green matrix are expressed through corresponding 
diagonal elements only by means of some linear operation. Specifically, from Eq. (0 the operator relation can be 
derived, 



G vl x = P„(l - R) 



(14) 



where both G VIX and G MAt are treated as matrices in extended x-coordinate space X. The operator R acts in the 
mixed coordinate-mode space M M constructed as a direct product of the space X and the truncated mode space which 
incorporates the whole set of mode indices except the unique index fx. In its turn, operator R is expressed as a product 
of operators and U, i. e. R = & V 'U, which are represented in M M by matrix elements 



{x,v\gW\x , ,v , )=G < ¥\x,x l )5 uu , 
{x,v\U\x' ,u') =U vv ,(x)8{x - x 1 ) . 



(15a) 
(15b) 



The function G^\x,x') in (|15a|l is the trial mode Green function which satisfies the equation resulting from 
provided that all inter-mode potentials are put identically equal to zero, 



dx 2 



+ k 2 u + iO-V u (x) 



GW(x,x') = 6(x-x') 



(16) 



The operator P„ in (|14fl is the projection operator whose action reduces to assigning the given value v (or fi) to 
the nearest mode index of arbitrary operator standing next to it (either to the left or right), without affecting the 
product in the X space. This operator may thus be thought of as contracting the action of scattering operator R from 
three-dimensional space M M to the unidirectional space X. 

Assuming \j! = \x and substituting into Eq. JSJ) all inter-mode propagators in the form (|14fl we arrive at a closed 
equation for the intra-mode (i. e. mode-diagonal) propagator G M/i (x,x'), namely 



^ + kl + iQ-V lx {x)-% 



G MM (x,a;') = 5(x - x') , 



for V/lx . 



Here. 



P„W(1-R)- 1 RP^ 



(17) 



(18) 



is the operator (in X space) potential, which in general case is highly non-local even though the magnetic-field-induced 
potentials may be completely non-existent. This potential bears a strong resemblance to the T- matrix which is well- 
known in quantum theory of scattering |34l l35l. Normally, this matrix is quite singular in the multi-channel case. 
However, in our approach, as was proven in |14| . the T-operator is properly regularized owing to separation of the 
intra- and inter- mode potentials in Eq. (|SJ) . We will refer below to this operator as the effective inter- mode potential. 

Bearing in mind the explicit form l|l(J|) of the potential V^(a;), in particular, the non-random nature of its "magnetic" 
part, it is worthwhile to renormalize unperturbed mode energy Q by introducing, in place of /c^, the magnetic-field- 
modified longitudinal energy 



2 ; 2 
>C. = K.. — 



H 2 



(19) 



As a result, in place of equations (|16fl and l|17fl we are led to analyze a couple of different, though equivalent, equations, 



dx 2 
dx 2 



+ x 2 v + iQ - V uv (x) 



G^\x,x')=5(x-x') , 



iO-V^ix)-^ 



G^(x, x 1 ) = 5{x - x 1 ) 



(20a) 
(20b) 



where solely the disorder potential V^^x) appears instead of the total potential (|10|) and the replacement has been 
performed k 2 ^ — > x^. 
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IV. SPECTRAL PECULIARITIES OF RENORMALIZED MODE STATES 



Before we proceed to the conductance calculation, let us examine fist spectral properties of the open quantum 
system in question, see Ref. [23 for more details. Although the true eigenstates are unknown in the presence of a 
random potential, equation l|20b|) permits of analyzing the particular mode state fi comprehensively. This can be 
done due to one-dimensional character of the equation of motion which can be regarded as completely decoupled from 
the remaining bath of mode states, whatever strength of the entangling potentials in Eq. (JSJ. 

Consider the operator in square brackets of Eq. Il20b|i . assuming equation H2()a[l to be ad interim solved (see 
Appendix) and the potential thus completely determined. Our aim is to examine the average value of this operator 
as it specifies the energy of the particular mode in the mean-field approximation. This approximation makes practical 
sense on condition that scattering produced by fluctuating parts of the potentials may be regarded as weak. As far 
as the disorder potential V^^x) is concerned, the weak scattering (WS) condition is normally cast to inequalities 



(21) 



where £ stands for the electron mean free path at a zero magnetic field. For the reference purpose, this path calculated 
from the model of the white-noise Gaussian-distributed potential, whose binary correlation function is (V(r)V(r 1 )') = 
Q5(r - r'), equals 4ir/Q. 

Given the external magnetic field, there appears an extra length parameter in the problem in question, specifically, 
the magnetic length, in terms of which it would be advisable to specify the condition where the magnetic-field- 
induced scattering could be regarded as weak. In our approach, this type of scattering is taken into consideration in 
two fundamentally different places. The first one is the "magnetic" part of the potential (|10|l . which is absorbed in 
the mode energy renormalization (Tffi . The other place is the mode- mixing potential (|18|l . In this study, we will focus 
on the limiting case of relatively weak magnetic fields, where maximal cyclotron radius of classical electron orbit, 
R c = kpl^, is large as compared to the quantum well thickness H. It was shown in [22j that under condition l|21(l 
of weak disorder-related scattering (WDS) in conjunction with the constraint 



(22) 



which will be hereinafter referred to as the condition for weak magnetic scattering (WMS ) , the norm of the inter- mode 
scattering operator R in l|18l) is small compared to unity, 



Rl 



1 



< 1 . 



(23) 



This limitation allows the inverse operator in l|18[) to be expanded in series and the operator potential 7^ to be 
simplified to the form 



f u » pM [v) up u 



(24) 



Unlike quasi-local intra-mode potential V^^x), the potential (|24|l possesses a nonzero mean value, even with 
no magnetic field. Bearing in mind the further applications of some perturbation theory, it is advisable to divide 
this operator into the mean and the fluctuating parts, = (7^) + A7^. In view of definition (|12|) . the average 

operator (7^) can be split, though quite conventionally, into the sum of quasi-local "disorder" and essntially non- 



local "magnetic" parts, viz. (T M ) = (Tj S) ) + (T^). The action of the operators (T^ y) ) and (f^ 1 ) on the Green 



(Q)\ 



-(B) 



function G fJ , fJ ,(x 1 x') is specified by equalities 



(t^ Q) )G^](x,x') 



E 



(25a) 



and 



[{t^)G^]{x,x>) = 



E 



2iH 



q( x ) C"H ;) 



(y) 



l 2 y dx 



H 2 
I 



( q( x ) 
2iH I ^ 



B 



sty d 



P y dX X 
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' B 



G^^(xi, x') 



(25b) 
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To proceed further with expressions Q25[l. it is necessary to specify the trial Green function G^\x, x') or, more 



precisely, its averaged value (G v [x,x')\. In doing so, we will use a somewhat simplified version of correlation 
relations for the disorder potential, namely 



(V(r)) = 0, 
(V(r)V(r')) = QW(x 



x')d(r x - r' ± ) 



(26a) 
(26b) 



As shown in the Appendix, under WS conditions which imply that WDS inequalities (|21|l and WMS inequality l|22|) 
are satisfied simultaneously, the averaging of the trial Green function results in the asymptotic expression 



(GiTW^^exp. 



1 



1 



\x — x 



(27) 



Here, L\0(y) are the forward (/) and the backward (b) scattering lengths of the trial mode v, which are determined 



from model i|26|) as 



45 
9Q 



{v) , ,_ 45 (2^)- 



(28) 



W(q) is the Fourier transform of the function W(x) from (|26b|) . Expression l|27|) is literally applicable to extended 
modes, i. e. to modes with kJ, > 0. with x 2 , < (evanescent modes) one should assume h v — i\kJ\ in Eq. (|27|l and 
both of the extinction lengths (|28fl to turn to infinity. 

With function l|27l) . the action of the operator (T^\ reduces to multiplying the operand in Eq. (|25a|l by the 



complex- valued disorder-related self-energy factor 0, , specifically, (T^\G^. 
4 S) = Ax*+i/Ti?\ with 



(x,x r ) 



^ 2) G MM (x,x'), where 



Q 



dq W{q + Xy) 
2ir q 2 — x 2 



» - 45 E — - **) + + *u) 



(29a) 
(29b) 



Symbol V in Eq. I|29a|) stands for the integral principal value, the bar over the summation index in (|29bfl signifies the 
summation over extended modes only, if any. 

The conditional character of the term "disorder self-energy" with reference to suggests that this factor is 

actually determined not solely in terms of the disorder potential, which results in the pre-factor of Q, but also is 
influenced by the magnetic field. The latter field renormalizes Wctvenunibers u and adjusts the number of extended 
modes, see below. Here it should be noted that the dephasing factor l|29b|) can only be nonzero if quench-disorder- 
induced scattering is allowed between different extended modes. Otherwise the coherence of the trial mode states is 
unaffected by static disorder. 

Apart from the disorder-induced fraction of the mode self-energy, the action of essentially non-local "magnetic" part 
of the operator results in additional self-energy term which will later on be referred to as magnetic self-energy, 
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(30) 



The principal difference between the disorder and the magnetic self-energies is that the former takes a non-zero value 
provided that the disorder strength (which is specified by the parameter Q) be finite. On the contrary, the magnetic- 
field-induced term only vanishes in the limit of the magnetic field equal to zero, remaining otherwise finite at 
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any disorder strength. By comparing imaginary parts of the terms 
of "magnetic" and "disorder" dephasing rates is evaluated as 



(Q) 



and sif^ one can determine that the ratio 
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(31) 



It is clear that under WMS condition (|22[) the magnetic-field-originated dephasing is negligible, whatever level of the 
disorder. One is thus led to conclude that strong inter-mode mixing resulting from the magnetic field only does not 
cause the mode energy levels to be considerably widened unless there exists a random potential due to some kind 
of disorder serving as a mediator for the magnetic-field-related dephasing effect. As far as mode entanglement is 
concerned, it is evident from i|29|) that the specific role of the magnetic field is to change the collective parameters 
of the electron spectrum, such as the mode content of the confined system of carriers and the mode density of states 
(MDOS). It is only in this indirect manner that the mode level width becomes dependent on the magnetic field. 

The number of extended modes in the quantum waveguide, which substantially governs dephasing rate i|29bfl . is 
specified by the requirement for mode energies to be positive- valued. Prior to analyzing in WS approximation the role 
of fluctuating potentials (x) and A 7^ one can notice that the energy of a specific mode depends not only on the 
confining potential but also on the magnetic field that likewise produces the confinement effect. As seen from <|19ll . 
the new "unperturbed" mode energy decreases steadily with the magnetic field growth, whatever the orientation 
within the plane of a 2D system. A decrease in that energy should result in truncation of the number of extended 
modes responsible for current transport and, therefore, in the conductance abruptly falling off. The latter effect is 
quite often interpreted in terms of electron effective mass enhancement [3(| . 

Although magnetic renormalization 119|) is of purely intra-mode nature, corrections of the same sign to the mode 
energy are produced by the inter-mode magnetic scattering as well, which is taken into account by self-energy 130fl . 
This is indeed the case as long as WMS condition holds true. Moreover, inter-mode magnetic scattering can lead to 
some anisotropy of the electron effective mass with respect to the magnetic field in-plane orientation, depending on 
symmetry properties of numerical coefficients in Eq. (|30|l . i. e. on geometrical symmetry of the quantum well. Note 
that the intra-mode magnetic term in (|19fl is entirely isotropic. 

As the magnetic field increases, the terms in square brackets of l|3Uf) . which are proportional to C^ v , will exceed 
the rest of the terms, thus resulting in the subsequent increase of the mode energy. This effect, however, corresponds 
to the range of relatively strong magnetic fields, where WMS inequality 12L'I) is violated. In this paper, this case will 
not be dealt with. 

In Figure [3 the numerical dependence of the extended mode number N c , commonly referred to as the number of 
effective conducting channels, is reproduced from Ref. |23|. The magnetic field assumed to be codirectional with the 



N„ 



25 



15 



kj.w/7^10.5 
— k F H/;r=5.l 



k H/tt=3 . 1 

F 

k H/tt=1 . 5 




10 



14 



FIG. 2: The dependence of the number of extended modes in the quantum waveguide on 
the inverse magnetic field [y is the Landau filling factor). 



axis of current flow is scaled as the Landau-level filling factor v = (fc^-Z^) 2 = kpR c - The collapse of the number of 
current-carrying modes with a growing magnetic field is apparent, regardless of the quantum well thickness H, the 
width W being held constant. Numerical considerations reveal that in-plane rotation of the magnetic field slightly 
changes the picture presented. This is consistent with the fact that the real part of self-energy l|30|l can reach, at most, 
the same value (on the order of magnitude) as the magnetic term in Ijl9(l , and is in agreement with experiments |37| . 
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We also reproduce from [27| the N c dependence on the effective thickness of the quantum waveguide, the latter 
adjusted in practice by depletion voltage. In the lowest magnetic fields (upper curve in Fig. the number of channels 
shows a near-linear increase with growing H. This corresponds to standard geometrical considerations applicable to 
systems of waveguide configuration. As the magnetic field is getting higher, the conventional geometric increase in 
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FIG. 3: The number of conducting channels vs the width of the near-surface potential 
well at different values of the Landau filling factor. The shaded area below the 
dashed curve covers the parameter region where the WMS criterion is violated. 



the number of channels slows down, gradually changing to a decrease in the conducting mode number. This unusual 
non-monotonicity of the mode content of the quantum system is due to essentially non-monotonic dependence on H 
of the magnetic-field-renormalized mode energy l|19l) . 

The parallel magnetic field effect on the number of subbands in a quantum well was previously noticed by some 
authors |38ll39|] with reference to semiconductor heterostructures. The reduction in the number of conducting channels 
with a growth in the magnetic field is certainly a quantum-mechanical effect pertinent to multi-electron systems. 
As a matter of fact, its nature is closely related to the Aharonov-Bohm (AB) phase incursion, which is known to 
effectively change carriers' energy HI |4l[ and has an impact on the interference patterns observed in experiments 
with AB rings. In the problem under study, unlike quasi-lD metal rings subjected to magnetic field, the AB phase is 
essentially different for electrons following different orbits. Therefore, it may seem to be difficult to take the magnetic 
field into account by conventional wave function phase renormalization. Yet this phase effect appears in the collective 
mode energies, affecting crucially the number of current-carrying modes in a laterally confined electron system. 

The apparent impact of the in-plane magnetic field on the number of extended modes inevitably occurs in the 
magnetic-field dependence of the mode dephasing rate. While in the limit of a large number of conducting channels 
the replacement of the sum in Eq. (|29b|l with the integral leads to the familiar quasi-classical formula, 



(32) 



where l/r^(0) = kpQ/^ stands for the dephasing rate due to disorder scattering in the absence of magnetic 
field 0|, in Fig. 0] we give the numerically obtained dependence on the magnetic field of two particular modes of the 
quantum waveguide under study. Well-known van Hove singularities of MDOS, which are directly related to step-wise 
change in the number of extended modes, are well pronounced on both of the curves. 

Besides, Fig.[S]illustrates the dephasing rate of the particular mode vs variable effective thickness of Q2D system of 
carriers. Availability of van Hove singularities similar to those depicted in Fig. ^compels one to conclude that by means 
of the orbital coupling to Q2D electrons the in-plane magnetic field has the effect which is, in a way, similar to that 
of the electrostatic confinement potential. At the same time, in contrast to the magnetic-field-controlled singularities 
shown in Fig. the oscillations of truly geometrical origin are noticeably more complicated. The distinction is 
attributed to a substantially different response of the effective mode energy (|T§|| to the magnetic field, on the one 
hand, and to size parameters of the confined electron system, on the other. However, it should be emphasized that in 
both of the graphs, Q] and [SI the reduction of the dephasing by the quenched disorder is clearly visible as the magnetic 
field grows. This fact is indicative of an increase in coherence of the electron transport in quench-disordered Q2D 
systems subjected to an external magnetic field. 
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FIG. 4: The dephasing rate (I29bll for two particular modes vs inverse magnetic field. 




FIG. 5: The dephasing rate vs the quantum waveguide thickness at different strengths 
of in-plane magnetic field. The broken fraction of the lower curve falls into the 
range of parameters where WMS condition 1122B is violated. 



V. MAGNETOCONDUCTANCE OF A Q2D SYSTEM 



By substituting the Green functions into Eq. (J2J) in the form of expansion in series over transverse Hamiltonian 
eigenfunctions we obtain the following mode representation of the conductance, 



g xx (L,B) 



dxdx' 
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(33) 



Here, 5 ViVj is the vector-argument Kronecker delta, z ViV . is the inter-mode matrix element of the z-coordinate operator 
which, given the model of the confining potential, assumes the value 
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Using interrelation (|14|l between off-diagonal and diagonal mode propagators, and estimate 1|23[) for the inter- mode 
scattering operator it is easy to make sure that the principal contribution to 133|l arises from the terms whose mode 
indices are coincident with one another. As a result, the conductance is expressed, asymptotically in WS limit, as 



9xx(L, B) 
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(35) 



12 



This formula, in conjunction with l|20|) . ultimately reduces the initially stated three-dimensional dynamic problem to 
a set of decoupled, strictly ID problems which in the general case are non-Hcrmitian. 

Under WS conditions, equation l|20b|l can be solved analytically with an arbitrary accuracy in fluctuating potentials. 
However, there is no need to overestimate the expected result. In Refs. [T4lll6|. it was proven that for inactive media 
with channel number N c > 1 the presence of dephasing term (|29b|) in the operator potential (|24|l enables one to 
neglect fluctuating potentials in Eq. (|20b|) . thereby seeking intra- mode propagators from the deterministic equation 



(^ + ^ + ! /^)g w (v') = ^~x') . (36) 

The solution to this equation, which obeys the openness conditions at the ends of the interval C = (-L/2, L/2), under 
WS conditions is written as 



G^.a/Jw^-exp^ix^-l/Z^la-a/l] , (37) 



where = 2xr At r^' ) is the length of the mode \i phase coherence. With function l|37|) substituted into Eq. I|35|l . 
the average magnetoconductance is given as 



<3 M (L,B))=£^ 



,(¥>) 

^ L 



1 /•* 

1 — ^— exp 




(38) 



Here, the terms corresponding to extended modes only are kept since the evanescent-mode Green functions, which 
are strongly localized in x-direction, are real- valued and cancel each other in Eq. (|35|) . Note that formula l|38l) can be 
viewed as describing the system conductance in the presence of both the random scatterers and the magnetic field, 
which enters implicitly through mode coherence lengths. 

From general result l|38|) . conventional limiting formulae for the conductance can readily be obtained. In particular, 
in ballistic limit L <C £ the dimensionless conductance becomes nearly equal to the number of open conducting 
channels, 

(g xx (L,B)) wJV c (B) . (39) 

This number, as it follows from the above analysis, is determined by both the geometric confinement of the electron 
system and the magnetic field, see Figs. |2 and |3] The perfect system conductance is thus expected to run in steps as 
a function of either depletion voltage or the value of the in-plane magnetic field. 

In diffusion limit L 3> £, if the potential well in the z-direction is wide enough to contain the large number of 
quantization levels in this direction, by replacing the sum in r.h.s. of (|38() with the integral we arrive at 

( 9xx (L,B)) » „ ( 9xx ( L ,0)) (l - ^ . (40) 

In the case of the zero magnetic field this result is coincident in form with classical Drude conductance If B ^ 0, 
the magnetoconductance is negative, being varied smoothly with the magnetic field, specifically, under the quadratic 
law. This is because the van Hove singularities of MDOS prove to be integrated out in such a rough calculation. 

However, these singularities are actually contained in the mode coherence lengths, as they are determined using 
the dephasing rate formula (|29b|) . They should appear both in Shubnikov-de Haas (i. e. magnetic-field-driven) 
oscillations of the conductance and in the conductance dependence on the quantum well width, which is normally 
tuned by depletion voltage. In Fig. [BJ the results numerically obtained from Eq. I|38() at several values of the diffusion 
parameter, L/£, are presented. MDOS singularities reveal themselves in the form of sharp dips placed close to those 
points where the number of conducting modes undergoes stepwise variations, i. e. to the thresholds of the transverse 
subbands, no matter what the disorder level may be. The disorder seems to manifest itself through the absolute value 
of the conductance. Note that in the case of relatively large disorder (large values of L/C) the conductance develops 
non-monotonically vs the magnetic field, even if MDOS singularities are smoothed out. When decreasing the mean 
free path, non-monotonicity becomes so apparent that it appears to be inadmissible to disregard this effect in the 
experimental data. In particular, with regard to this analysis it would be tempting to revise the observed positive 
magnetoresistance which is frequently attributed to spin properties of 2D systems |42l 1 13) . 

To make a comparison with the magnetic-field run, in Fig. the conductance is shown versus the width of the 
potential well forming a quasi-2D quantum waveguide. Here, MDOS singularities are equally more pronounced, which 
demonstrates the change in the number of conducting channels. Meanwhile, in the latter case these singularities are 
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completely anticipated from the very outset, since the number of channels is normally associated with size quantization. 
The peculiar feature should be noted in Fig. as against the dependence on the magnetic field. The curves in Fig. 
tend to become non-monotonic, on average, as the magnetic field grows. Evidently, this non-monotonicity accounts 
for the non-monotonic dependence on H of the mode eigen-energy <|19fl . 
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FIG. 7: The conductance vs the width of the near-surface potential well at different values of the in-plane magnetic field. 

In panel (b), the region selected by the rectangle in panel (a), with only two curves left, is zoomed in to make the 
dips arisen due to MDOS singularities more clearly visible. Similarly to Fig. |3J the shading below the dashed curve 
signifies the region where criterion 1221 is violated. 

It is instructive to dwell upon the physical nature of dips in Figs. Eland [7\ All of them are positioned in the vicinity 
of the points corresponding to opening/closing of the conducting channels. For one thing, as the electron waveguide 
gets thin or the magnetic field grows, the closing of the channel must result in a bench-like fall of the conductance 
since each of the conducting channels is expected to bring in exactly one conductance quantum. At the critical point, 
the edge extended mode is transformed to the evanescent one, which is localized at a scale of the mode wavelength 
and in this way carry no current in an infinitely long system. Indeed, this picture is demonstrated by the upper, 
"ballistic" , curve in Fig. [5] 

For another, in approaching the transformation point (subband threshold) the MDOS of the edge mode diverges 
whereas its mode velocity tends to zero. The capacious and slow edge mode serves as a destructive sink for dynamic 
electrons, leading to a decrease in the conductance. It is clear that the dips can arise only when we deal with an 
imperfect system of carriers, where scattering is allowed from all remaining extended modes to the slow critical mode. 
The conductance in the bottom of the dip may thus reach, hyperbolically, a nearly zero value, which can hardly be 
grasped in real experiments because of various unaccounted extra factors. 
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VI. DISCUSSION AND CONCLUDING REMARKS 

We have demonstrated that the appreciable localizing effect produced by relatively weak in-plane magnetic field on 
2D electron and hole systems can be reasonably interpreted in the context of Fermi liquid theory applied to spinless 
electrons which reside in an open near-surface potential well of finite rather than zero width. For a relatively weak 
magnetic field, its coupling to the carrier orbital degree of freedom, even though it is rather weak from semi-classical 
point of view, proves to have substantial influence on the carrier spectrum and, hence, on the conductance. 

The conclusion about strong sensitivity of the carrier spectrum to the in-plane magnetic field is made from the 
analogy of tightly gated solid-state systems to classical wave-guiding systems of planar, though three-dimensional, 
configuration. The mode content of these well-known objects is quite sensitive to the anisotropy in their cross-section, 
that is to the applied gate voltage as far as electron devices are concerned. A remarkable feature of the latter type 
of systems is that in case the magnetic field is applied to the finite length, these systems can be thought of as being 
subjected to both the disorder potential, whose correlation length can be arbitrary, and the additive deterministic 
strongly non-local "magnetic" potential barrier. Scattering parameters of this barrier are specified by the magnetic 
field strength and orientation, on the one hand, and by the length of magnetically biased section of the quantum 
waveguide, on the other. The effect produced by this barrier results from the mismatching of electron spectra in 
the inner and outer parts of the quantum well, the inner part representing the finite-length electron system under 
consideration. 

If the electron waveguide cross-section is strongly anisotropic, the bulk of transverse modes in the electron spectrum 
proves to be efficiently transformed from extended to evanescent type as the magnetic field grows slightly. This is 
because electron scattering from side boundaries of the confining potential well is, in the strict sense, specular if the 
boundaries are considered as completely inhibiting the transverse current flow. In view of this fact, it seems to be 
erroneous to assess the magnetic field effect on the electron system by simply estimating the variation of the electron 
trajectory portions between successive collisions with quantum well boundaries. It can be easily verified from Eqs.© 
and 119|) that the mode truncation effect of the magnetic field is the more significant the larger is the aspect ratio of the 
waveguide cross-section, given the dimension H . For H being small enough, such that only modes with quantization 
parameter m = 1 in the corresponding direction can be regarded as extended, the total number of extended modes in 
the quantum waveguide is mainly determined by the larger cross-section dimension, W in the case under study. Owing 
to this, even slight alteration of the magnetic field can transform a considerable number of modes from extended to 
evanescent type, thus leading to significant reduction of the conductance, even though it might have a large (ballistic) 
value in the absence of the magnetic field. 

The mode truncation effect of in-plane magnetic field appears to be quite similar to that of truly geometric con- 
finement of the electron system. The closing of each of the current-carrying modes, which makes itself evident in the 
form of conductance jumps by exactly one conductance quantum in a perfect system at a zero temperature, should 
be regarded as a true quantum phase transition |44j . This statement is substantiated by the indisputable fact that 
there exists a well-defined correlation length in the vicinity of a closing point, whose role is played by the wave length 
of the marginal extended mode, which tends to diverge as the critical point is approached. The closing of the last 
conducting mode by means of the in-plane magnetic field may thus be regarded as the magnetic-field-driven MIT. 

It should be noted that in the proximity to the MIT the conductance is not quite precisely described by the 
present theory, since it is hard to satisfy WMS conditions over the corresponding range of system parameters. At the 
same time, closer examination of magnetic-field-originated items in Eq. (|20b|) reveals that the above-described mode 
truncation effect and, hence, the very fact of the existence of magnetic-field-driven MIT, is robust. 

In this paper we have not focused our attention upon the conductance anisotropy with respect to in-plane orientation 
of the magnetic field. The corresponding analysis is straightforward. It consists in comparing the real part of 
anisotropic self-energy (|30|l and the magnetic part of the trial mode energy (|19|l . Under WMS condition (1221 , both of 
these energy items are of the same order of magnitude, save the case of the magnetic field oriented nearly parallel to 
the direction of current. We thus expect the conductance anisotropy to be expressed, at most, moderately. Moreover, 
in view of numerical coefficients in Eq. (|30[) being model-specific, the conductance anisotropy in practice should be 
governed substantially by the confinement potential profile. It should be noted that assuming the magnetic field to 
grow beyond WMS limit, the isotropic terms proportional to l^ 4 can be made dominating in self-energy (|3L)fl . thereby 
removing the in-plane anisotropy of the conductance. This range of the magnetic field strength needs to be given 
a special consideration because the "magnetic scattering" , which is an important aspect of our approach, cannot be 
regarded as weak in this particular case. 

Yet another relevant remark should be made concerning the model of the confinement potential adopted in this 
study. In some papers where laterally confined electron systems arc dealt with (sec, c. g., Ref. 25]), this potential is 
taken as a quadratic function. Thus modelled confinement seems to be beneficial from the technical point of view, as 
it enables one to account for the magnetic field non-perturbatively, the corresponding transverse eigen-states being 
known as Fock-Darwin levels [3^, |3t| . It is natural that the precisely zero width of those levels in the absence of any 
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disorder implies the entire lack of dephasing due to the magnetic field only. In this case, the magnetoconductance 
should not exhibit a dip structure because the latter results from MDOS singularities arising exclusively in the presence 
of the disorder. 

In the domain of weak magnetic fields (in a sense of inequality l[22J) the quadratic confinement can hardly be 
substantiated. Therefore, the issue of the magnetic-field-induced mode entanglement and the related dephasing of the 
mode states might seem to be quite topical in this limiting case. Note, however, that although rectangular confinement 
possesses appreciable inter-mode magnetic scattering, it proves not to result in widening the transverse quantization 
levels unless some kind of disorder is also taken into account. The exact form of transverse eigen- functions is of 
no fundamental significance for the development of transport theories in mode representation, but the mere fact of 
transverse energy quantization does matter. This prompts us to expect that the magnetic field alone cannot give 
rise to noticeable decoherence of electron states for any hard- wall model of the electron confinement. The specific 
form of the confinement potential can only alter the arrangement of transverse quantization levels and, consequently, 
influence the coherence properties of the electron system indirectly. Of primary value is the dephasing resulting from 
scattering caused by some kind of random (i. e., in a sense, uncontrolled) potential, no matter static or variable it 
may be. 
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APPENDIX: TRIAL MODE GREEN FUNCTION IN THE PRESENCE OF MAGNETIC FIELD 
1. Reduction of the boundary- value problem to auxiliary Cauchi problems 



Prior to obtaining the proper boundary conditions for Eq. 121) aj , which would conform to the physical definition 
of the system openness at the end points x — ±L/2, note that the problem governed by this equation is absolutely 
identical to that of point-source-radiated classical waves in an ID random medium. Bearing this in mind, we will 
extend the driving equation from the finite interval C to the whole x-axis, writing it down in symbolic form 



2 



A 

— + k 2 + iO - V(x) 
ax* 



G^(x,x') = S(x - x') , (A.l) 



where the potential V(x) is assumed to be the finite-support function which consists of two fundamentally different 
items. The regular component of this potential is V^ re9 '(x) = y^r (l — ^ m -i ) 0(L/2 — \x\), with 6{x) being the 

Heaviside unit-step function, whereas the additional random disorder term is V i -- dls \x) = Vun(x)d(L/2 — \x\). For the 
sake of clarity we will omit mode index fi, implying that both the wave number k in (|A.1|) and n below are related 
exactly to this mode. 

Since there is a need to perform configurational averaging over the potential U( rf ") (x) it is worthwhile to express 
the solution to Eq. (|A.lfl in terms of wave functions of causal type rather than functions that comply with the initially 
stated boundary-value (BV) problem. This is achieved by employing the formula 

G( v Hx,x') = ^[i J+ (x)^(x')9(x-x / ) + iP + (x')iP-(x)6(x' -x)] , (A.2) 

where ip± (x) are two different solutions of homogeneous equation i|A.l|) with the boundary conditions specified for each 
of them at only one end of the coordinate axis, viz. x — > ±oo, depending on the sign index, W is the Wronskian of the 
above solutions. With this representation, the trial propagator itself satisfies, as it must, the initial BV problem and 
the averaging procedure can be accomplished mathematically correctly, with due account of multiple backscattering. 
Inasmuch as the support of the potential in Eq. I|A.1|) is bounded, functions ip±{x) may be sought in the form 

ip ± (x) = TT ± (x)e l ^ ±x - L ^ - i 1± {x)e- l ^ ±x - L / 2) at \x\ < L/2 , (A.3a) 

ip±(x) = c±e *H±x-L/2) at > L j 2 ; (A. 3b) 

where x 2 = k 2 — V^ reg \ Under WDS condition l|21|) . envelope functions n±(x) and J±(x) in l|A.3a(l can be regarded 
as smooth factors as compared with near-standing fast exponentials, which leads to the following coupled dynamic 
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equations, 



±tt±(x) + ir](x)Tr±(x) + (l(x)~f±(x) = , 
±j'±(%) - iv{x)l±(x) + C±(x)ir±(x) = . 



(A.4a) 
(A.4b) 

Random functions rj(x) and (±{x) in (|A.4I) are constructed as narrow packets of spatial harmonics of the disorder 
potential, 



C±(x) = — / -V^(t)eM^(±x-L/2)} 
ZKJ x -l 21 



(A.5a) 
(A.5b) 



where spatial averaging is over the interval 21 of an arbitrary length intermediate between "small" length scales x^ 1 
and r c , on the one hand, and the large scattering length (to be determined self-consistently) , on the other. These 
limitations ensure smoothed random "potentials" rj{x) and C±( x ) to provide the harmonics ±>r forward and backward 
scattering, respectively. 

By joining the solutions (|A.3a|) and (|A.3b|l at the end points of interval C we arrive at the necessary BC for envelopes 
n±(x) and r y±(x), viz. 



7r±(±L/2) =const , 
7±(±L/2) =& B \ ± (±L/2) . 



The quantity 



.fc- 



(A.6a) 
(A.6b) 



(A.7) 



as it follows from (|A.3a|l . is the amplitude reflection coefficient from the interface between the magnetically biased and 
unbiased regions of the extended ID quantum waveguide. This reflection will be hereinafter referred to as "magnetic 
scattering" associated with the above introduced potential V^ re9 \x). It can be easily verified that under constraint 
(|22|1 the reflection parameter i|A.7|) is modulo small as compared to unity. 

By substituting functions ip±( x ) i n the form l|A.3af) into i|A.2fl . the trial Green function whose both coordinate 
arguments are located inside magnetically biased interval C is given as 



G (y ) (x, x') = £i (x, x'y^ x -^ + g 2 {x, x')e-' l ^ x - x "> + G 3 {x, x')e l ^ x+x "> + £ 4 (z, z'K"*^* 

where the last two terms in r.h.s. have emerged as a direct consequence of the mode backscattering. In 
envelope functions are 
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where the notations are used 



A(x) = [l + r + {x)r_{x)e 2l * L ] 
T±(x) = j±(x)/tt±(x) , 
6± = 0[±(x-x')\ . 



(A.8) 
smooth 

(A.9a) 
(A.9b) 
(A.9c) 
(A.9d) 



(A.lOa) 
(A.lOb) 
(A.lOc) 



The functions T±(x) play a particular role in the averaging technique. It can be deduced from Eq. (|A.3a(l that these 
functions are nothing but the reflection coefficients of spatial harmonics ±x which are incident at the point x onto 
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the disordered layers whose end coordinates are x and ±L/2, respectively. These reflection factors are well-known to 
obey the Riccati-type dynamic equations |45|. 



±^j^ = 2i V (x)T ± (x) - ( ± (x) + C±(x)r%(x) , 
which in our case must be supplied with "initial" conditions (one-side BC) stemming from ()A.6|I . viz. 

r±(=LL/2) = . 



(A.ll) 



(A.12) 



2. Configuration averaging of the function 

The averaging technique as applied to functionals of smoothed potentials (|A.5|) was elaborated in Refs. 0,EHE3|- 
It relies basically on the fact that under WS conditions the functional arguments (|A.5|) can be regarded as Gaussian- 
distributed random variables 48] . Taking account of correlation equalities i|2ti|) , the binary correlators of these effective 
potentials, which would suffice to be allowed for, can be cast to the form 



(i]{x)ri(x')) 



-Ft(x- x') , 



(<; ± (x)a(x')) = --F l (x-x') 



(A.13a) 
(A.13b) 



where Lf and Lb are the forward and backward scattering lengths given in (|28|) for the particular mode v. The 
function 



Fi(x) 



2tt {ql) 2 21 



21 



9(21- \x\) 



(A.14) 



is sharp on a scale of significant change in smooth amplitudes of electron wave functions, and therefore it can act as 
the under-limiting (5-function when averaging amplitude factors Gi(x, x') in IjA.Sj) . Prior to averaging these factors, it 
is advantageous to absorb, wherever possible, the forward-scattering field rj(x) by going over from functions r y±{x), 
tt±(x) and C±i x ) to phase-renormalized functions 



7± ( x ) = 7± ( x ) ex P 
tt± ( x ) = 7T ± ( x ) exp 
C±(x) = (±(x)exp 



±L/2 

±i / dxir)(xi) 

±L/2 

! ' / dxi77(xi) 

±L/2 

±2i / dxirj(xi) 



(A.15a) 
(A.15b) 
(A.15c) 



An important thing is that correlation relation (|A. 13b|) remains unchanged after renormalization (|A.15c|l . 

In view of short-range correlation of random functions l|A.5p and due to the causal nature of functionals being 
averaged, the averaging of functionals with different sign indices in (|A.9|1 can be done separately. By applying the 
Furutsu-Novikov formula for Gaussian random processes |45| , we first average the equation for renormalized reflection 
factors r±(x) = 7 Y±(x)/tt±(x), i. e. 



± 



dr ± (x) 

dx 



-C±(x) + Ci(x)r 2 ± (x) 



With condition l|A.12Jl taken into account this yields 



(r±(a;)) =ft (B) exp 



1 (L 



Tx 



(A.16) 



(A.17) 



In view of WMS-motivated smallness of the reflection parameter H^ B \ this implies that only the terms that do not 
contain the factors of T± (x) and 7± (x) can be kept in i|A.9|) when averaging Green function 
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Next, in order to average the ratio 7r±(x')/7r±(x) consider its Fourier transform over x' coordinate, 

r±L/2 



& ± \x,q)=±[ dxx^^exp 
Jx 7r±(a0 



-iq(x — X\) + ik\x — X\ | ± i / dx2?7(x2) 



(A.18) 



where we have separated the forward-scattering random field tj(x) in the form of a phase factor. Averaging over this 
field yields readily 



exp 



±i / dx 2 r](x2) 



\x - Xi\ 
„ 6XP ^ 2L f 



Then, after averaging over f](x), the function (|A. 1 8|) comes to obey the equation 

ix T iq) <$ (±) (a:, q)) v ~C± ± (x)($ (±) (x, q)) n 



d(& ± \x,q)) n _ j /j. 
dx 



2L f 



(A.19) 



(A.20) 



which is to be solved jointly with Eq. I|A.16[) . Averaging of Eq. HA.20|) over the effective back-scattering field (±(x) 
results in the dynamic equation 



d($W(x,g)) _ i 



dx 



(i (±) (^,9)) 



(A.21) 



which is to be solved with obvious "initial" condition, (<&( :iz >(±L/2, q)) — 0. The solution to Eq. (|A.21|I yields 
1/1 1 



2 \Lf Lb 
wherefrom the result is obtained 



ik =F iq 



1 — exp 



- I — — 

2 VL/ + L b 



ik =p zg 



(i/2Tx) 



(A.22) 



(Gi(x,x'))e i ^ x - X H (G 2 (x,x')) 



'\\ P - i>c ( x - x ') 



e v ' — — exp 
2k 



\k\x - x'\ - ^ + ) ./■ 



(A.23) 



The terms (jA.9cj) and l|A.9d(l can be averaged in the same manner. The result is much more cumbersome, namely 



(Gz{x,x')) = — - — exp 
2k 



ixL — L [ 1 

Lf 2Lb 



0(x — x') exp 



[x — x) exp 



xfl 1 \ x[_ /_3 1 
2 \T S + T b ) ~ ~2 \ T } + T b 

x' ( 1 1\ x / 3 1 

^ \T f + Tb) ~ 2\T S + T b 



(A.24a) 



(£ 4 (x,x')) 



2^ 



exp 



" ! l T JT, 



(x — x ) exp 



(x — x) exp 



2 \L 





X 


/ 3 


1 


— ) 


+ 2 




'T 


M 


x' 


(3 _ 


1 


W 


+ 2 


U/ 


h L 6 



(A.24b) 



Fortunately, in the WMS limit, because of the factor being small, these terms may be omitted, leading finally 
to result l|27}. 
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